THE INTEGRALS IN GRADSHTEYN AND RYZHIK. 
PART 6: THE BETA FUNCTION. 



VICTOR H. MOLL 



Abstract. We present a systematic derivation of some of the definite integrals 
in the classical table of Gradshteyn and Ryzhik that can be reduced to the 
beta function. 



1. Introduction 

The table of integrals [2] contains some evaluations that can be derived by ele- 
mentary means from the beta function^ defined by 

(1.1) B{a,b)= [ x^-^l-xf-^dx. 

Jo 

The convergence of the integral in ()1.1|) requires a, b > 0. This definition appears 
as 3.191.3 in 0. 

Our goal is to present in a systematic manner, the evaluations appearing in 
the classical table of Gradshteyn and Ryzhik [2 , that involve this function. In 
this part, we restrict to algebraic integrands leaving the trigonometric forms for a 
future publication. This paper complements [3j that dealt with the gamma function 
defined by 

/•oc 

(1.2) r(a) := / x''-^e-''dx. 

Jo 

These functions are related by the functional equation 

n -^^ w w r(a)r(&) 

r(a + b) 

A proof of this identity can be found in [1] . 
The special values r(n) = {n — 1)1 and 

for n E N, will be used to simplify the values of the integrals presented here. Proofs 
of these formulas can be found in [3] as well as in Proposition 12.11 below. 
The other property that will be employed frequently is 

(1.5) r(a)r(l-a) = 

sm na 
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The reader will find in a proof based on the product representation of these 
functions. A challenging problem is to produce a proof that only employs changes 
of variables. 

The table [2j contains some direct values: 

Jo {l-x)P ~ sinpTT 
is 3.192.1 and is evaluated by identifying it as B{p+ 1, 1 —p). Formula 3.192.2 is 

dx IT 



(I 7) 

^ ^ Jo {l- x)P+^ sinpTT 

has the value B{p + 1, -p) = T{p + l)r(-p). Next, 3.192.3 is 

(1.8) /'^^ ' 



xP+^ sinpTT 
and the change of variables t — 1/x in 3.192.4 produces 

(1.9) r{x - l)P-i/2^ ^ r <-P-i/2(i _ 
Jl X Jo 

and this is 

(1.10) B(i_p,i+p) =r(i-p)r(i+p) = , 

\- - / \z / \z / COS p-TT 

as stated in [2]. 

Let 6 = i in to obtain 



_ n _ r(a) 



(1.11) / ^^^B{a,^) = 

The special 
respectively. 



lo ^ r(a+i)' 

The special values a — n + \ and a — n + ]^ appear as 3.226.1 and 3.226.2, 



2. Elementary properties 

Many of the properties of the beta function can be established by simple changes 
of variables. For example, letting y = 1 — x in (|l.ip yields the symmetry 

(2.1) B{a,b) = B{b,a). 

It should not be surprising that a clever change of variables might lead to a 
beautiful result. This is illustrated following Serret 4j. Start with 

B{a,a) = [ {x - x^)"-'^ dx 
Ja 

- 2/ dx. 

Jo '- ^ 

The natural change of variables v — ^ — x yields 

/•1/2 

(2.2) B{a,a)^2 / (l - v^)""^ dv. 

Jo 

The next step is now clear: let s — Av'^ to produce 

(2.3) B{a,a)^2^-^''B{a,^) . 
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The functional equation (ll.3|l converts this identity into Legendre's original form; 
Proposition 2.1. The gamma function satisfies 

. ^ r(2a)r(i) 

In particular, for a = n G N, this yields (|1.4p . 

3. Elementary changes of variables 

The integral (jl.ip defining the beta function can be transformed by changes of 
variables. For example, the new variable x — t/u, reduces (II. ip to 



(3.1) / e-'{u-ty-' dt = u''+''-'B{a,h), 
Jo 

that appears as 3.191.1 in [2^. The effect of this change of variables is to express 
the beta function as an integral over a finite interval. Observe that the integrand 
vanishes at both end points. Similarly, the change t = {v — u)x + u maps the interval 
[0, 1] to [u, v]. It yields 

(3.2) / {t-u)''-\v-t)''-Ut^ {v-u)''+'"'^B{a,b). 

J u 

This is 3.196.3 in The special case w = 0, w = n and a = h = n~\-\ appears 
as 3.193 in [2] as 

Several integrals in can be obtained by a small variation of the definition. For 
example, the integral 



(3.4) 



/ {l-x'^f-^dx = -B{l/a,b) 
Jo a 



can be obtained by the change of variables t = x"". This appears as 3.249.7 in [2] 
and illustrates the fact that it not necessary for the integrand to vanish at both end 
points. The special case a ^ 2 appears as 3.249.5: 

(3.5) /V - x^)'-^ dx = \B (i, fe) = 22^-25(&, 6), 
Jo 

where the second identity follows from Legendre's duplication formula 
The change of variables t — cx produces a scaled version: 

(3.6) / (c" ~t'')''-Ut=-c''^''-^'>+^B(l/a,b). 
Jo a 

The special case a — 2 yields 

rc 26-1 

(3.7) / (c2-t2)''-idt = ——5(1/2,6). 



2 

The choice b ~ n + ^ appears as 3.249.2 in [5]: 
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Similarly 3.251.1 in ^ is 
(3.9) / x^-^l-xy-'dx^-Bl-.b 



Jo " 
The change of variables t ~ 1/x converts (|l.ip into 

/oc 
t-''-''(t-l)''-idt==B(a,6). 

Letting t — xP yields 



(3.11) 



l<00 -I 

The special case ly — b and /i = p{l — a — b) is 3.251.3: 

poo -, 

(3.12) / a;''-i(2;P-l)"^' dx= -B(l-i.-/i/p,i/), 

Ji P 



4. Integrals over a half-line 

The beta function can also be expressed as an integral over a half-line. The 
change of variables t — x/{l — x) maps [0, 1] onto [0, oo) and it produces from p.ip 

(4.1) B{a,b) 



(1 + *)"+^' 

In particular, if a -|- 6 = 1, using (|1.3p and (jl.Sp . we obtain 

°° dt _ TT 

1 + i sin 7ra 
This can be scaled to produce, for a > and c > 0, 



(4.2) 







(4.3) / —^^^d'-^ foroO 







X + c sm na 

that appears as 3.222.2 in [5]. In the case c < we have a singular integral. Define 
b——c>0 and s = so now we have to evaluate 

poo a-l J 

(4.4) / = -b'^-' / -. 

Jo 1 - s 

The integral is considered as a Cauchy principal value 



(4.5) / = lim , , , , 

Let y = 1/s in the second integral and evaluate them in terms of the beta function 
to produce 

1 / r(a) r(l-a-e) 



(4.6) /=limer(e)x , „, . , 

e-^o e \I(a-fe) 1 (1 — a) 

Use L'Hopital's rule to evaluate and obtain 

r'(a) r'(i - a) 



(4.7) I = - 



r(a) r(a) 
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Using the relation r(a)r(l — a) = Trcosec 7ra, this reduces to TrcotTra. Therefore 
we have 

(4.8) rx-^Ux^ —{-cY~^ forc<0 

Jo X -\- c tanTra 

The change of variables x = e~* produces, for c < 0, 

(4.9) / = -ncot{t,n) {-cr-\ 
J-oo e * + c 

The special case c = —1 appears as 3.313.1: 

/■°° e"^* dt 

(4.10) / = 7rcot(/i7r). 

J-oo 1 - e ' 

We now consider several examples in [2] that are direct consequences of (14. 3p and 
(j4.8p . In the first example, we combine (|4.3p with the partial fraction decomposition 

1 _ 1 / 1 1 

{x + a)(x + b) h — a\x + a x + b 

leads to 3.223.1: 

(4.12) / , , = Tr^K"' - 6^-')cosec(^A^). 
7o (a; + fo)(a; + a) ^ ^ ^ 

Similarly, 

(4.13) ^ ' « + ^ 



a; + 5 a; — a (a — x) (b + x) 
leads to 3.223.2: 

f°° X^~^ dx TT 

(4.14) + - ^ (^^-^ -sec(^.) + a-^ cot(M.)) 

using gJl) and gH). The result 3.223.3: 

/■°° x^'-^dx , ,a^-i-6^-i 

(4.15) / — = 7rcot(^7r)- 



foUows from 
(4.16) 

Finally, 3.224: 
(4.17) 

follows from 
(4.18) 



(a — x)(6 — x) h — a 

1 1/1 1 



(a — x)(6 — x) a — h\h — x a — x^ 
(x + b)x^~^dx TT fa^b c — b 



(x + a) (x + c) sin(/x7r) \ a — c 

x + 5 6 — a 1 b — c 1 



(x + a)(x + c) 



c—ax+a c—ax+c 



We can now transform (|4.ip to the interval [0,1] by splitting [0,oo) as [0,1] 
followed by [1, oo). In the second integral, we let t = 1/s. The final result is 



"1 



(4.19) B{a,b)= / — '-—j-dt. 

^ ' Jo (1 + *)"+'' 
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This formula, that appears as 3.216.1, makes it apparent that the beta function is 
symmetric: B{a, b) — B{b, a). The change of variables s ~ 1/t converts (|4.19|) into 
3.216.2: 

(4.20) B{a,b)= / ———-ds. 

It is easy to introduce a parameter: let c > and consider the change of variables 
t — cx in (|4.ip to obtain 

(4-21) / -^^-^ = c-'^B{a,b), 

Jo [l + cx)''+'> 

that appears as 3.194.3. We can now shift the lower limit of integration via t — x+u 
to produce 

(4.22) / {t-uy\t + v)-''-''dt = {u + v)-''B{a,b), 

J u 

where v — 1/c — u. This is 3.196.2, where v is denoted by /?. Now let 6 = c — a in 
the special case u = to obtain 



(4.23) / {t-u)''-H-''dt:^u''-''B{a,c-a). 

J u 

This appears as 3.191.2. 

We now write (|4.ip using the change of variables t = x'^. It produces 

x'"'"^ dx 1 



Jq (1 + .1'^)"+'' C 

The special case c — 2 and a = 1 + /i/2, 6=1 — /i/2 produces 3.251.6 in the form 

X'^+^ dx ^ A^TT 

^ ' ' Jo (1 + 2:2)2 4smfnT/2- 

Now let & = 1 — a and choose a = p/c to obtain 

(4.26) / ——^ = -By-, = - cosec(7rp/c). 

Jo l + x'^ c \c c J c 

This appears as 3.241.2 in [2]. 

Similar arguments establish 3.196.4: 

(^•2^) / ia-bx)ix-ir =-b'°''''^^^^ [b-a 

Indeed, the change of variables t = x — \ yields 

^ ' ' Ji [a-bx){x-ir ~ Jo [{a-b)-bt] t-' 

and scaling via the new variable z ^ bt/{b — a) gives 

. . „„x . dx 1 f b [°° dz 

(4.29) 



{a-bx){x-lY b\b-aj Jq (l + z) z"^ 
The result follows from (|4.ip and the value 

(4.30) B{iy, l-iy)= r(i^)r(l - i^) = —~ 



Sm TTl' 
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The same argument gives 3.196.5: 



(4.31) / — — — = — cosec(i^7r) 



_^ (a ~ bx){l — b \a — b 



5. Some direct evaluations 



There are many more integrals in P] that can be evaluated in terms of the beta 
function. For example, 3.221.1 states that 

(5.1) r = - hf-^ cosec Trp. 



X - 



To establish these identities, we assume that a > 6 to avoid the singularities. The 
change of variables t = {x — a)/{a — b) yields 

(5.2) r^^^^^^ia-bY-^ r^-^^^ 



X — b Jq 1 + 1 ' 



and this integral appears in 
Similarly, 3.221.2 states that 

(5.3) / ; = -n[b — ay cosec np. 

J^oo x-b 

This is evaluated by the change of variables y — [a ^ x)/{b ~ a). 

The table contains several evaluations that are elementary corollaries of (|4.ip . 
Starting with 

(5-4) / , „ =Bia + l,b-a-l) = 



we find the case a — p and = 3 in 3.225.3: 



(5.5) 



/■°° X- 

Jo (T" 



Pdx _T{p+l)T{2-p) _p{l-p) TT 
+ a;)3 ~ f(3) ^ 2 sin(p7r) ' 



using elementary properties of the gamma function. 

The change of variables t — 1 + x converts (|5.4p into 
(5.6) r = i?(a + M - a - 1) - + ^) - « - ^) 



The special case a = p — 1 and b — 2 gives 

(5.7) r ^'"Vr'"^' = r(p)r(2 - p) = (1 - p)r(p)r(i - p) - " 



sin(p7r) 

This appears as 3.225.1. Similarly, the case a = l—p and 6 = 3 produces 3.225.2: 



(5.8) 



^^^^ = ^(^-^)\(^+^) = 1,(1 -P)r(p)r(i -p) = 

i3 r(3) 2^^ ^' ^ ^' 2sin(p7r) 
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6. Introducing parameters 

It is often convenient to introduce free parameters in a definite integral. Starting 
with (|4.ip . the change of variables t = ^x'^ yields 



This formula appears as 3.241.4 in with the parameters 
(6.2) a=— ,b = n + l — — ,0^ = 1^, u = q, and v = p, 

V V 

in the form 

x^-^ dx _ 1 (vY'" r(^/t/)r(n + i- iijv) 

/o {p + qx"Y^^ ~ \q) r(n+l) 

This is a messy notation and it leaves the wrong impression that n should be an 
integer. 



I' 

Jo 



• The special case v = c = 1 and b = p + 1 — a produces 
/ , f'-^dt 1 , 

Jo (1 + ut)P+^ u"- 

This appears as 3.194.4 in P], except that it is written in terms of binomial coef- 
ficients as 

/ , f°° t°-~^dt , TT /a - 1\ 

We prefer the notation in 



The special case v = c = \ and b = 2 — a produces 



Using (|1.3p and (|1.5p yields the form 

'"^ f'-^dt _ (1 -a)7r 
(1 + ut)^ u°- sinTTo 



(6.6) 

This appears as 3.194.6 in [2] 



• The special case u = v = 1 and c = q, and choosing a = p/q and b = 2 — p/q 
yields 3.241.5 in the form 

xP-'^ dx ^ q-p TT 
^ ' ' Jo (l + a;9)2 ^2 sin(7rp/g)' 

• The special case c = 1 and a — m + 1, b — n — m — ^ produces 
(6.8) / - = -B{m + l,n-m- ^) 

Using (|1.3p and (|1.4p this reduces to 

(aq] r ^"''^^ _ m!n!(2n-2m-2)! ^,„,^, ^"^-"+i/2 

^ ■ ^ 7o („ + „f)«+5 (n-m-l)!(2n)! ' 
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for TO, n e N, with n> m. This appears as 3.194.7 in [2]. 



(6.10) / = -S (a, i - a) 



• The special case u = v = \ and 6 = i — a yields 

'■"^ r^-^dt _ 1 
/o \/TTF " c 
Writing a — p/c we recover 3.248.1: 

fOC ,p-l 7, -, 

• Now replace v by in (|6.ip . Then, with u = 1, a = i, c = 2, so that ac = 1 
and 6 = n — i we obtain 

_ 1 

/o (t;2 + f2)rx - 2^-1 

This can be written as 

ffin^ /°° _ - 1/2) 

^ ' Jo (t>2 + i2)n 2r(n)t;2»-i 

that appears as 3.249.1 in [2]. 



(6.12) / 7ZT^ = ^B{ln-l] 



The special case w 1, c = 2 and & = ^ — a in (|6.ip yields 

(l + ut2)«/2 ~ 2^ 



foo ,2a-l Jf 1 

(6-14) / ^ - TT^^ («. T - a) 



Now a = 1 /2 gives 



(6.15) 



It is curious that the table [2] contains 3.249.8 as the special case m l/(n — 1) of 
this evaluation. 

• We now put u = v = 1 and c = 2 in ()6.ip . Then, with b = 1 — v — a and 
a = /i/2, we obtain 3.251.2: 

^ ' Jo (l + i2)i-^ 2 V2' 2 

• We now consider the case c = 2 in (|6.ip : 

^ ' Jo (w + Ut2)a+b ^ ^ 

The special case a — ni + ^ and b — n ^ m + ^ yields 

t^"^dt r(TO+ l/2)r(n- m+ 1/2) 



(6-18) („^^^2)«+l 2u"+l/2„n-m+l/2r(n+l) ' 

and using (|1.4p we obtain 3.251.4: 

t2"rfi _ 7r(2TO)!(2n-2TO)! 



(6-19) (w + ut2)"+i 22"+iTO!(n-TO)!n!M'"+i/2x,"-™+i/2' 

for n, TO G N with n > m. 
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On the other hand, if we choose a = m + 1 and b = n — m we obtain 3.251.5: 
^2m+i^^ _ r(m + l)r(n-m) _m!(n-m-l)! 



(6.20) 



(w + ut2)n+i 2w'"+iw"-™r(n + 1) 2n!-it'"+it;'' 
Several evaluation in [2J come from the form 

1 



(6.21) 



6-1 



dt 



q 



-B{a,b), 



obtained from (II. ip by the change of variables x = f^. 

• The choice a = I + p/q and = I — p/q produces 

(6.22) / tP+'^-\l-t'>)-P/Ut^-B(l + -,l- 
Jo Q \ 1 



pn 



cosec 



Ptt\ 
1 J 



This appears as 3.251.8. 

• The choice a — 1/p and b = 1 — 1/p gives 
(6.23) 

This appears as 3.251.9. 



/' Wf-i(l - x«)-i/P dx^-B(-.l--] ^ - cosec ( -] 
Jo 1 \P' Pj Q \PJ 



• The reader can now check that the choice a = p/q and b = 1 — p/q yields the 
evaluation 



Jo q \q qj q 



(6.24) 

This appears as 3.251.10. 

• Putting V = I and = — a in (|6.ip we get 



cosec 



piT 

q 



(6.25) 



oo .ac— 1 



dt 



, = B(a, V — a). 

(1 + cu"^ ' 



Now let a — r/cto obtain 



(6.26) 



oo ,] — 1 



dt 



(1 + uV'Y cu 



1 /? 

.r/c V ( 



This appears as 3.251.11. 

• We now choose b ~ 1 — 1/q 'm (|6.2ip to obtain 

6.27 / ^-B\a,l-- 

^ ' Jo </T^ q \ q 

Finally, writing a = c — (m — l)/q gives the form 



(6.28) , ^ 

Jo q 

The special case q — 2 produces 



/' 

Jo 



1 



= -B I C+ ,1- - 

q q q 



(6.29) 



71372 2^1^+2 2.2. 



r(c- 



2r(c + i-f) 
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In particular, if c = n + 1 and ?7i = 1 we obtain 3.248.2: 
(6.30) 



/o 2r(n + 3/2) (2n+l)! 

Similarly, c ^ n and m — yield 3.248.3: 

t^"dt _ TT (2n)! _ TT f2n 

In the case g = 3 we get 

6.32) / ^-B{r+ ,1-- 

^ ' Jo 3^33' 3 

This includes 3.267.1 and 3.267.2 in [2J: 



3V3r(i)r(n+l) 



3^ 

^t^'^-^dt (n-l)!r(|) 



^JT^ 3r(n+|) 
The latest edition of 2 has added our suggestion 

as 3.267.3. 



</T^ 3r(n+i) 



7. The exponential scale 

We now present examples of (jl.H) written in terms of the exponential function. 
The change of variables x = e~'^* in (|l.ip yields 



(7.1) / e-''\l-e-''y-' dt= -B[-,b 



(7-2) / , ...,,a+, =-B{a,h). 



This appears as 3.312.1 in :2]. On the other hand, if we let x = e in (|4.ip we 

get 

_ 1 

(l + e-^*)° ^ 

This appears as 3.313.2 in [2j. The reader can now use the techniques described 
above to verify 

r°° e~^^^ dx \ { vw 

(7-3) (^bATj— =«exp S(aAi,z.-a/i), 

that appears as 3.314. The choice 6 = 0, = 1 and relabelling parameters by 
a—Xjq and \x = v yields 3.311.3: 

(7.4) / _ = -B -, 1 - - = -cosec — 

using the identity B{x, 1 x) — 7rcosec(7ra;) in the last step. This is the form given 
in the table. 
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The integral 3.311.9: 

(7.5) / J — — nb^^^ cosec(/^7r) 

J^oo o + e ^ 

can be evaluated via the change of variables t = and ()4.2p to produce 



(7.6) / = 6^ 







1 + t 

8. Some logarithmic examples 

The beta function appears in the evaluation of definite integrals involving loga- 
rithms. For example, 4.273 states that 



(8.1) 



The evaluation is simple: the change of variables x — ut produces, with c = v /u, 

(8.2) / = / \uP-^t{\nc~\ntY~^ —. 

Ji i 

The change of variables z — give the result. 
A second example is 4.275.1: 



(8.3) 



/' [(- Inx)"-! - x^-\l - x)^-i] = [Tip + q)- T{p)] , 



r(p + g) 
that should be written as 

(8.4) I [{-\nxY~^ -xP-^{l~xy-^]dx^T{q)-B{p,q). 
Jo 

The evaluation is elementary, using Euler form of the gamma function 

(8.5) r{q) = f (-lnx)«-^ dx. 

Jo 

9. Examples with a fake parameter 

The evaluation 3.217: 

, , f°° f bPxP-^ (l + bx)P~^\ , 

9.1 / -; — dx = 7rcot7rp 

^ ' Jo \{l + bx)P bP-^xP J 

has the obvious parameter b. We say that this is a fake parameter in the sense that 
a simple scaling shows that the integral is independent of it. Indeed, the change 
of variables t = bx shows this independence. Therefore the evaluation amounts to 
showing that 

(9.2) r ( - ii±^)!li^ dt^TT cot TTp. 

To achieve this, we let y = 1 /t in the second integral to produce 

<-°° tP-^-" dt f-^dt 



^'■'^ ^^oj, il + t)p J, (1 + t)-.- 
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The integrals above evaluate to B{p — e, e) — B{e, 1 — p~ e). Using 
(9.4) B{a, b) = Ij"'^^^'!] and r(a)r(l - a) - ^ 



r(a + b) sin(7ra) 

this reduces to 

^ ' ' \ er{p)r{p + e)sm{TT{p + e)) 

Now recall that 

(9.6) limer(e) = l 

and reduce the previous limit to 

^^■^^ ^ ^ T^ip) sin(7rp) H \ ^^^P " '^^^P + '"'(^(P + ~ ^'(^^ 
Using L'Hopital's rule we find that / = 7rcot(7rp) as required. 

The example 3.218 

9.8) / ^ dx^TTCOtTip 

Jo ia + x)PxP 

also shows a fake parameter. The change of variable x — at reduces the inte; 
above to 

(9.9) / , — ^ — dt^TT cot np 

^ ' Jo {l + t)PtP 

This can be written as 
The result now follows from (19.21). 



10. Another type of logarithmic integral 
Entry 4.251.1 is 

f°° Ina; , i^b"'^ „ , 

(10.1) / ; dx — (In 6 — TT cot TTfl) . 

Jo X + b sm 7ra 

To check this evaluation we first scale hy x — bt and obtain 

10.2) / dx = b°-^Hnb hfe^^W dt. 

Jo x + b Jo ^ + t Jo 1 + t 

The first integral is simply 

r°° t^-i ^ 

(10.3) / ——=B{a,l^a)=T{a)T{l-a) = - . 

Jo 1 + r sm TTo 

The second one is evaluated as 

/■°° t"^^ Ini , 9 cosTTo 

10.4) / dt^-n^ — ^ 

^ ^ Jo l + t sin2(7ra) 

by differentiating (|4.ip with respect to a. The evaluation follows from here. 
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11. A HYPERBOLIC LOOKING INTEGRAL 

The evaluation of 3.457.3: 
(11-1) r . . = f ^, a Ina, 

is done as follows: write 

, „x T 1 X dx 
11-2) I=— 1 \ ^ 

and let t = ae^ to produce 

The change of variables s — yields 

_ 1 f"" Insrfs Ina r°° s^'/^-^ ds 

^ ~4^Ja (1 + 2^ Jo {1 + sY ' 

The first integral vanishes. This follows directly from the change s i-^ 1/s. The 
second integral is the beta value indicated in the formula. 
In particular, the value a — 1 yields 

r°° xdx _^ 
7_oo cosh'' X 
Differentiating with respect to /i produces 

oo 



(11.6) / zlncoshxdx = 0, 

J —CO 

that appears as 4.321.1 in [5]. 
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